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Abstract 

We study the statistical behavior of two out of equilibrium systems. 
The first one is a quasi one-dimensional gas with two species of parti- 
cles under the action of an external field which drives each species in 
opposite directions. The second one is a one-dimensional spin system 
with nearest neighbor interactions also under the influence of an exter- 
nal driving force. Both systems show a dynamical scaling with domain 
formation. The statistical behavior of these domains is compared with 
models based on the coalescing random walk and the interacting ran- 
dom walk. We find that the scaling domain size distribution of the 
gas and the spin systems is well fitted by the Wigner surmise, which 
lead us to explore a possible connection between these systems and 
the circular orthogonal ensemble of random matrices. However, the 
study of the correlation function of the domain edges, show that the 
statistical behavior of the domains in both gas and spin systems, is not 
completely well described by circular orthogonal ensemble, nor it is by 
other models proposed such as the coalescing random walk and the 
interacting random walk. Nevertheless, we find that a simple model 
of independent intervals describe more closely the statistical behavior 
of the domains formed in these systems. 
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1 Introduction 



In this article we study of the statistical behavior of two non-equilibrium 
systems. The first one is a quasi one- dimensional gas introduced in Ref. pp. 
There, the authors studied the biased diffusion of two species in a fully pe- 
riodic 2 x L rectangular lattice half filled with two equal number of two 
types of particles (labeled by their charge + or — ). An infinite external field 
drives the two species in opposite directions along the x axis (long axis). The 
only interaction between particles is an excluded volume constraint, i.e., each 
lattice site can be occupied at most by only one particle. 

The second system is a one-dimensional spin system introduced in Ref. [2] , 
where the authors consider a chain of L Ising spins with nearest neighbor 
ferromagnetic interaction J. The chain is subject to spin-exchange dynamics 
with a driving force E that favors motion of up spins to the right over motion 
to the left. 

Both systems evolve with formation of domains which grow in time. For 
intermediate times where the size of the domains is much smaller than the 
total size L of the system, the domain size distribution exhibit a dynamic 
scaling. In this work we will be interested in the statistical properties of 
these domains in the scaling regime. 

For each system there are two kind of domains. For the gas system, there 
are domains filled with particles and empty domains. For the spin system 
there are domains of up spins and domains of down spins. We are interested 
in the statistical properties of these domains regardless of their type (filled 
or empty, up or down). Let us define P^ n \S,t) as the probability density 
function that the distance between the external borders of n + 1 consecutive 
domains is S at time t. Let (S) be the average of S and the relative spacing 
between borders is s = S/ (S). The scaling probability density function is 
defined as 

p^(s) = (S)P^(s(S),t), (1) 

in particular p(°>(s) is the scaling nearest neighbor edges distribution, i.e. the 
domain size distribution. In the scaling regime, the scaling probability den- 
sity functions p( n >(s) do not depend on the time t. 

Numerical simulation of both systems performed in [TJ [2], [3] show that 
their scaling domain size distribution function is well fitted by 

/»( s ) = | se — 2 /4. (2) 



2 



This distribution is known in the random matrices theory as the Wigner 
surmise and it describes the spacing distribution between eigenvalues in the 
Gaussian and circular orthogonal ensembles (GOE and COE respectively) [I]. 
This fact suggest a possible connection between the statistical behavior of the 
borders of domains in non-equilibrium systems and the eigenvalues in random 
matrices in a similar way as it occurs with the eigenvalues of the Gaussian 
orthogonal/unitary ensemble and the vicious random walk [21 El E]. 

The main objective of this paper is to explore this possible connection 
between non-equilibrium domain systems and random matrices, in particular 
the circular orthogonal ensemble. In sections [2] and [31 we summarize the main 
results found in Ref. pQ on the quasi one- dimensional gas, and in Ref. [2J [5] on 
the spin system. In section m we recall some facts about the random matrix 
theory, relevant for our purposes. In section [51 we compare the statistical 
behavior of the domains for each system with the statistical behavior of the 
eigenvalues of random matrices from the circular orthogonal ensemble. Also, 
we will consider a simple model of independent intervals, which turns out 
to better describe the statistics of the domains in both the gas and the spin 
systems. 

2 Quasi one-dimensional gas 

This system was described briefly in the introduction, more details are found 
in Ref. pQ . The system evolves in time according to the following dynamical 
rules: 

1. L particles are randomly inserted in a 2 x L rectangular lattice, ^ par- 
ticles (+) and ^ particles (— ), the remaining sites are empty. Periodic 
boundary conditions are imposed in both directions of the lattice. Let 
the x axis be the long axis of length L. 

2. Two neighbor sites are chosen at random. The contents of the sites 
are exchanged with probability 1 if the neighbors sites are particle- 
hole, but if they are particle-particle the content are exchanged with 
probability 7. The exchanges which result in +/— particles moving in 
the positive/negative x direction are forbidden due to the action of the 
external field. 

3. A time unit correspond to 2L attempts of exchange. 
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2.1 Qualitative results 



The system was studied by computer simulation in Refs. [U [8]. We also 
performed several simulations of this system with a parameter 7 = 0.1. The 
following qualitative results are found from the simulations: 

• For low 7 values the system remains homogeneous, i.e., the system 
evolves without domain formation. 

• For high 7 values the system evolves as follows. For early times small 
domains form everywhere due to mutual obstruction of the opposite 
species. After a long time, the system settles into a non-equilibrium 
steady state (NESS) in which only one macroscopic domain survives, 
so that, its length is about -| and it contains almost no holes. In the 
low density region the traveler particles leak out of the domain at one 
end and later rejoin it at the other end. The charge distribution of this 
macroscopic domain is not trivial. 

• For intermediate times, when the average size of the clusters is much 
smaller than the size L of the system, the system shows dynamical 
scaling. 

2.2 Quantitative results 

Following Ref. we can map the quasi one- dimensional lattice into a one- 
dimensional lattice. This approximation described in Ref. [Tj is a coarse 
grained description. For any configuration on the 2 x L lattice we construct 
an effective one-dimensional one, with occupation numbers zero or one on 
a L sites line, as follows. At each site i, we assign if there are 5 or less 
particles in the 10 sites around it, including the z-th column of the original 
lattice. We assign 1 otherwise. In the coarse grained configuration a domain 
is a simple consecutive sequence of ones and its size is just the length of this 
string. 

Under this coarse grained description, the authors of Ref. [1] considered 
the "residence distribution" of the clusters P(S,t), defined as P(S,t) = 
SP(°'(S,t) in terms of the domain size distribution p(°'(S,t). Using the 
residence distribution, they found [H [8] that the average length of domains 
S = J2s SP(S, t) I P(S, t) grows in time with an exponent of at least 0.6. 
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The numerical results from the simulations of Ref. |IJ, show that the scaling 
residence distribution p(s) = SP(sS, t) is well fitted by 

p(s) = - 2 x\-^'- . (3) 

Remembering that p(s) is, by definition, proportional to sp ( -°\s), this implies 
that the scaling domain size distribution p(°\s) is properly fitted by the 
Wigner surmise for the orthogonal ensembles of random matrices, Eq. ([2]). 
We performed our own simulations to confirm this, the results will be shown 
in section [51 

In Ref. P, the authors proposed an approximate model to describe this 
system, in the coarse grained description, using the coalescing random walk 
(CRW) where the particles execute independent random walks, suffering a 
fusion reaction (A + A — > A) when two particles meet. This approximation 
is useful because the full hierarchy of correlations for this diffusion-limited 
reaction has been solved, see Ref. [91 [TOj [HI [12] . This approximation is done 
under the following assumptions. 

• The numbers of domains is equal to the number of walkers in the CRW. 

• The walkers describe random walks over a one dimensional ring with N 
sites, in that way, N its equal to the sum of the length of the domains. 

• The distances between walkers are equal to the length of the domains. 

• When two walkers meet they coalesce, representing the disappearance 
of a domain. 

• The traveling particles of the original model are neglected in the CRW 
description. 

Under this approximation the spacing distribution is the same found as in 
the coarse grained description, but the average length of domains grow in 
time with an erroneous exponent of 0.5. 

In order to solve this problem, in Ref. [T], the authors introduced the 
model of the interacting random walkers (IRW), where neighboring walkers 
experience stronger attraction with decreasing separation. The probability 
of a walker moving to the right is 

— (-ra 2 ) 
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(4) 



where l r is the distance to its right nearest neighbor and q s is the probability 
of staying. In the same way, the probability of moving to the left is 



where // is the distance to its left nearest neighbor and q s = 1 — q r — qi. 
The factor 0.06 (S) is chosen in Ref. [1] in order to give the best fit to the 
simulation data and the dependence with l~ 2 is justified in Ref. pQ. In this 
model, S grows with an exponent of 0.6, see Ref. pp. The spacing distribution 
is still given by 

We performed a few simulations of this interacting random walk changing 
the constant 0.06 in Eqs. (jl]) and (J5]). This has the effect of changing the 
growth exponent of the domains, however the scaling domain size distribution 
p(°)(s) remains unchanged, it is still given by the Wigner surmise fl2]). 

3 One-dimensional spin system 

In the introduction we briefly presented this system, for more details see 
Ref. 013]. The lattice has a length N with Nfi spins up ("+") and N(l-fi) 
spins down ( "— " ) with < fi < 1. Periodic boundary conditions are imposed. 
The spin-flip events are: 

1. + + <->+-+- A = AJ-E. 

2. --++«->- + -+ A = 47 + E. 

3. + + -+ «-> + - ++ A = —E. 

4. - + <_>__+_ A = —E. 

where the transition probability rate for a process from left to right is propor- 
tional to | (l — tanh (^)). The constant J is the nearest neighbor coupling 
between spins, E is the energy associated to an external field which drives 
the up ("+") spins to the right and the down ("— ") spins to the left, and T 
the thermal energy (temperature times Boltzmann constant). 

As explained in [2], the microscopic dynamics of the lattice of spins may 
be mapped onto one for an array domain dynamics, which provides a good 
approximation in the regime T -C E <C J ■ In this approximation, domains of 
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up spins move spontaneously to the right, and domains of down spins move 
to the left. The rates for such processes are independent of the domain size. 
The algorithm used for the numerical simulation is the following: 

1. Set up a random array of alternating down and up spins, with /iiV spins 
up and (1 — n) N spins down. 

2. Choose a domain at random. 

3. If the domain is down, move it to the left (i.e., reduce the size of its 
left neighbor by one and increase its right neighbor by one), otherwise 
move it to the right. 

4. If one of the neighbor domains is zero size, then the domain is removed. 

5. Update the clock by 1/number of domains. 

6. Repeat steps 2-6. 

This simplified system allows much better statistics than it would be 
possible with the true microscopic system [2] , and its simulation is more easy 
to implement and it reproduces the same asymptotic behavior. 

3.1 Qualitative results 

From the simulation, the following qualitative results are observed: 

• For early times, little domains form everywhere. 

• At intermediate times, some domains disappear while other domains 
grow. In this time regime, the system shows dynamical scaling. 

• For later times, only two macroscopic domains remain, which move in 
opposite directions. The system falls in a non-equilibrium steady state. 

3.2 Quantitative results 

In Ref. [2], the authors derive an analytical solution in the case fi — > 0, 
nevertheless we can derive another exact solution for the dynamical domain 
model with [i = 0.5. Let P^(n,m,t) be the frequency of finding a domain 
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of spin up with length n and another domain of spin down with length m at 
time t. The master equation for p(°\n,m,t) is 



dPW(n,m,t) = p(0)( n + ljmjt) + p (0) (n _ ljmjt )_ 2p (0)( njmjt ) 

(JjL 

+ P (0) (n,m + l,t) + P (0) (n, m - 1, t) - 2P (0) (n, m, t) . 

(6) 

Let Ax be the distance between nearest neighbor sites of the lattice. In the 
continuous limit, Ax — > 0, the above probability is a function of S = nAx, 
R = mAx. It satisfies the master equation 

dP®(S,R,t) _ d 2 P^(S,R,t) d 2 P^(S,R,t) 

dt ~ OS 2 + OR 2 U 

If at t — all domains have length So, the above equation is subject to 
the initial condition p(°\S, R, 0) = 5(S — Sq)5(R — So) and the boundary- 
conditions P^°^(0,P, t) = and p(°\S, 0, t) = 0. The appropriate scaling 
solution for this differential equation is 

2 2 

p ( j (s,r) = — sre 4 e 4 , (8) 

as a consequence, the length of domains of spin up and of domains of spin 
down have the same probability distribution function found in the quasi one- 
dimensional gas. The above equation can be written as the product of two 
distributions with the form of the Wigner surmise (121) . 

The fact that the domain size distribution p(°>(s) for a domain of up or 
down spins is well described by the Wigner surmise © was also confirmed 
in Ref. [21 [3] from the numerical simulation results. 



4 Random matrices 

Historically, the initial motivation to introduce random matrices was to study 
the statistics of the energy levels of some quantum systems, in nuclear physics 
initially, then for quantum systems when their classical counterpart is chaotic, 
for a review see Ref. [3]. Some local properties, such as the spacing between 
consecutive high energy levels are well described by the random matrix the- 
ory. For quantum chaotic systems with time-reversal invariance and integral 
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total angular momentum, the high energy levels, for a given value of the 
spin and the parity, have spacing distributions and correlations that are well 
described by the Gaussian orthogonal ensemble (GOE) of random matrices, 
when the energy levels are rescaled such that the average interspacing is 
unity. For large matrices, the (properly rescaled) statistics of the eigenvalues 
of the GOE in the bulk is the same as the one for the circular orthogonal 
ensemble (COE). In this work we preferred to work with the circular orthogo- 
nal ensemble rather than the Gaussian orthogonal ensemble since the former 
one is periodic and that way we avoid boundary condition problems (in the 
GOE the statistics of the eigenvalues near the edge are different than in the 
bulk). 

The COE is an statistical ensemble of N x N unitary symmetric matrices 
S invariant under the transformations S — > W T SW where W is any N x N 
unitary matrix. The eigenvalues of S are of the form e k , with k — 1, . . . , N. 
The probability density function for the eigenvalues is 

P N ^9 1 ,...,9 N ) = ^- J] \e»>-e*'\ P (9) 



N/3 



l<k<j<N 



with (3 = 1. The partition function is Z Np = (2n) N T(l + /3N/2)/(T(l + 
j3/2)) N [I]. It should be noted that the eigenvalue probability density func- 
tion of the COE is equal to the Boltzmann factor of an equilibrium system 
of particles moving on a circle and interacting through a log-potential, the 
two-dimensional Coulomb potential: 

P Np (e 1 ,...,e N ) = ^-e-^'-^ (io) 

with 

V N (6 1 ,...,6 N )= -ln|e^-e^|. (11) 

l<k<j<N 

Up to an additive constant, the log-potential is scale invariant, thus in 
these log-gases the density plays only the trivial role of fixing the length scale 
of the problem. In the other systems we study here, it is important to rescale 
the lengths in order to have the interparticle average spacing equal to one. 
For the non-equilibrium systems, this has been possible due to the fact that 
they exhibit a scaling regime. 

The COE is a homogeneous system, the density of eigenvalues is constant 
on the circle, equal to 2ir/N. It is convenient to rescale the lengths such 
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that the average spacing between eigenvalues is unity. Let Xk = N&k/(2ir) 
be the rescaled eigenvalues (note a small abuse of language here, since the 
eigenvalues are really e tdk not 9k)- Analytic expressions for the eigenvalue 
spacing distributions are known, however they are complicated, see Ref. [I]. 
It turns out that a very accurate approximate expression for the nearest 
neighbor spacing distribution p(°\s) is given by the Wigner surmise ([2]). In 
the limit of large matrices, N — > oo, the correlation function between two 
(rescaled) eigenvalues x and x' is given by 



with r = \x — x'\ and s(r) = sin(7rr)/(7rr). 

5 Results and discussion 

5.1 The nearest neighbor spacing distribution and the 
Wigner surmise 

We performed simulations for the gas, the spin, the CRW, the IRW and the 
circular orthogonal ensemble (COE) of random matrices with the following 
parameters. For the gas, the spin and the CRW systems, we used a lattice 
with 1000 sites, and for the IRW the lattice had 500 sites. The random 
matrix simulations were performed with 200x200 matrices. For the non- 
equilibrium systems, the simulations were carried out long enough to reach 
the scaling regime. The initial density of particles in the CRW and the IRW 
was 1/2 and 2/3 respectively. In the scaling regime, we build histograms for 
the domain sizes (i.e. the nearest neighbor spacings for the CRW, IRW and 
COE). The data to build the histograms was taken at three different times 
for each system in order to verify the existence of a proper scaling regime. 
These times were T = 1000, T = 1500 and T = 2000 for the gas system; 
T = 10, T = 18 and T = 34 for the spin system (the time unit of each system 
was defined in the previous sections); T = 50, T = 100 and T = 200 for the 
CRW, and T = 100, T = 150 and T = 200 for the IRW. For the IRW and 
the CRW systems, the Monte Carlo time unit corresponds to a trial to move 
all the particles of the system. To have appropriate statistics, we performed 
20000 simulations (realizations) of the gas, the spin and the COE systems, 
and 50000 realizations of the CRW and the IRW systems. 




(12) 
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Figure [T] shows the scaling domain size distribution p(°\s) for the gas and 
the spin systems, obtained from the simulations. It is compared with the 
scaling nearest neighbor distribution of random walkers in the IRW and the 
CRW, and the scaled nearest neighbor distribution of eigenvalues of the cir- 
cular orthogonal ensemble, which we also obtained from simulation of these 
systems. We also present the Wigner surmise, Eq. (j2j), which is know to 
numerically reproduce accurately the nearest neighbor distribution of eigen- 
values of the circular orthogonal ensemble (it is exact for 2x2 matrices). 

As we can see in figure [TJ the Wigner surmise reproduces correctly the 
nearest neighbor distribution p^ (s) for the gas, the spin, the IRW, the 
CRW and the COE systems. This suggest a possible connection between the 
random matrices and the out of equilibrium systems. 

Motivated by this possible link between the non-equilibrium statistical 
systems and the random matrix theory (i.e. the equilibrium log-gas system), 
we explore the possibility if not only the nearest neighbor distribution is the 
same for all these systems, but also the other spacing distributions p^(s) 
for n > and the pair correlation function g(r). The approach here is to 
consider the edges of the domains in the gas and spin system as fundamental 
entities, and to study their spacing distributions and correlation function, 
in the scaling regime. With this study we wish to know if the similitude 
between these systems goes beyond the nearest neighbor distribution, or if it 
is simple a coincidence that the nearest neighbor distribution is the same for 
all systems. This also can give an indication about how much information 
on each system is contained in the nearest neighbor distribution. 

From our simulations, we computed the other spacing distributions p^ (s) 
and the correlation function g(r) of all these systems in order to compare 
them. The results are presented in the following subsections. 

5.2 The quasi one-dimensional gas and its approximate 
models: the IRW and the CRW 

In Ref. [T] , the domain size distribution of the quasi one-dimensional gas was 
described by the approximate models of random walkers, the IRW and CRW. 
We now compare the other spacing distribution functions and the correlation 
function. The results are shown in figure [2J The CRW and IRW reproduce 
correctly the nearest neighbor distribution p^ (s). But, both the IRW and 
the CRW fail describing the other spacing distribution functions p^ 1 (s) with 
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Figure 1: Nearest neighbor spacing distributions for the gas, spin, IRW, CRW 
and COE and Wigner surmise. 

n > 0, and, as consequence, they also fail describing the two-point correlation 
function g (r), since the correlation function can be expressed as 



We notice in particular the fact that the correlation function of the gas system 
shows an oscillation around r = 1 (and probably more for higher values of 
r). On the contrary, the CRW and IRW models show no oscillation in the 
correlation function g{r). For these systems g(r) seems to be a monotonic 
function. 

The spacing distributions p( n >(s) for n > 1 of the IRW and CRW are 
somehow similar, but show higher maxima values than those of the gas sys- 
tem. Furthermore, the maxima of the IRW and the CRW are located at 
higher values of s than the corresponding ones of the gas system. 

We conclude that the CRW and IRW systems are not statistically equiv- 
alent to the quasi one-dimensional gas, although these models were proposed 
in Ref. pQ to describe (approximately) the quasi one- dimensional gas. 
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Figure 2: Spacing distributions and the pair correlation function in the gas 
system, IRW and CRW for a lattice with N = 500 and N = 1000 sites 
respectively. 

5.3 The gas and spin systems and the circular orthog- 
onal random matrix ensemble (COE) 

Now, in figure El we compare the spacing distributions p ( - n \s) and the corre- 
lation function g(r) of both the gas and spin systems with the ones for the 
random matrix COE. 

The correlation functions g(r) of the gas and the spin system seem to be 
very similar, almost identical, with only small differences of a few percents in 
relative difference. They both exhibit the oscillation near r = 1 mentioned 
above. On the contrary, the correlation function for the COE does not show 
any oscillation, and it is very different from the correlation function of the 
gas and spin system. 

The spacing distributions p( n >(s) are somehow similar for the gas and the 
spin system, with some small differences between them. But, in any case, 
they differ much from the ones of the COE. The distributions p^ n \s), for 
n > 1, are more "disperse" in the gas and spin systems than the circular 
orthogonal ensemble, i.e. they have a width at half height larger than the 
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ones for the COE. Also, the maxima for each distribution of the gas and spin 
system are located at smaller values of s than for the COE. The maxima 
values are also smaller for the gas and spin system than for the COE. 

It is possible to conclude that the gas and spin systems have a similar 
statistical behavior but this behavior is different from the COE. They only 
coincide in the nearest neighbor distribution. 
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Figure 3: Comparison between the gas, spin system (/x = 0.5) and COE for a 
lattice with N = 1000 sites. In the COE we took 20000 matrices of 200 x 200. 

In the previous subsection, we have shown that the CRW and the IRW 
systems do not describe properly the statistics of the spacing distributions of 
the gas system. But do they have a similar behavior as the one of the COE? 
In figure HJ we compare the correlation functions and spacing distribution 
functions of both the CRW and IRW systems with COE. 

The correlation function g(r) of the CRW and IRW differs from the one 
of the COE. Although the three correlation functions are monotonous (they 
have no oscillations), around 1 < r < 2, the correlation function of the COE 
is smaller than the one of the IRW and the CRW. The first two distribution 
functions p^(s) and p^(s) of the three systems, IRW, CRW and COE, are 
very similar. However from n > 2, the COE spacing distributions p( n \s) 
start to differ from the corresponding ones of the CRW and the IRW. The 
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spacing distributions of the CRW and the IRW seem more "disperse" (as 
defined above) than the ones of the COE. The maxima values for the spacing 
distributions in the COE are larger than for the IRW and CRW. Also, these 
maxima are located at smaller values of s for the IRW and CRW than for 
the COE. 

We conclude that the CRW and the IRW systems are not statistically 
equivalent to the COE. On the other hand the CRW and the IRW show 
very similar spacing distribution functions and correlation functions. Thus, 
another important conclusion is that the interaction proposed in Ref. [I] for 
the IRW do not change considerably the statistical behavior of the domains 
of the system from the CRW in the scaling regime. The interaction between 
neighbors, Eqs. (j3J and change the growth exponent of the domains. 
But, once rescaled, the statistics of the walker spacing distributions and 
correlations of the CRW and IRW are very similar. 
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Figure 4: Comparison between the CRW, IRW and COE for a lattice with 
N = 1000 sites and in COE we take 20000 matrices of 200 x 200. 
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5.4 The independent interval approximation (HA) 

The results presented above show that the statistics of the edges of the do- 
mains of the out of equilibrium gas and spin systems, properly rescaled in 
the scaling regime, are not described by the COE, thus are not described by 
an equilibrium statistical system of particles with a log-potential interaction. 
Furthermore, for the gas systems, the approximate models of CRW and IRW 
do not reproduce correctly higher order spacing distributions and the cor- 
relation function. In this section, we propose a model that describe better 
the statistics of the domains of the out of equilibrium systems in the scaling 
regime when properly rescaled. 

We have found numerical evidence that suggest that the statistical behav- 
ior the domains of both the gas and spin systems can be fairly well described 
by the independent interval approximation (HA), where the correlation be- 
tween domains is neglected. With this simplification, it is possible calculate 
the distributions j/ n )(s), for any n > 1, see Refs. (131 [HI [15l [161 EH]- For 
example, for n — 1, the probability to have two domain edges separated by 
a distance s, knowing that there is one edge in between at a distance r from 
the first edge is p^(r)p^(s — r), if the domains are considered independent. 
Therefore, 

pO-\s)= f drp {0 \r)p(°\s-r). (14) 



It is useful to introduce here the Laplace transform of the distribution func- 
tions 

poo 

pW(Z) = / p^\ s )e- ls ds. (15) 



Taking the Laplace transform of equation (I14p . we have 

p (1) (0 = (p W) (0) a - (16) 

If the Wigner surmise ([2]) is assumed for p(°\s), its Laplace transform is 

= l-le l2/ *eric(^J , (17) 

where erfc(z) = (2/y/ir) e~* 2 is the complementary Gaussian error func- 
tion. Using the inverse Laplace transform of expression (|T6l) or from a direct 
calculation of Eq. (fT4"|) . we find 



P W (s) = ^e^ Us + V2e^ (-4 + tt, 2 ) erf ( ^ ) ) , (18) 
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with erf(z) = 1 — erfc(z). 

More generally, under this approximation, the Laplace transform of p^ (s) 
is simply 

p<">(0 = (^ O) (0) n = fl-^erfcf-^V . (19) 
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Since the pair correlation function can be obtained from Eq. (1131) . then g(r) 
is given by the sum over convolutions of p^ (r) , see Ref . [13] . In terms of the 
Laplace transform g of the correlation function, we have 

5 (i) = g r(i)) * = _^W_ j (20) 

with p^°'(l) given by Eq. (TlTl) . We have 

1 - le l2 / w erfc ( 4=) 

m = — — #■ (21) 

Ze' 2 /- erfc f^J 

The inverse Laplace transform of the above expression can be computed 
numerically to obtain g(r). 

In this independent interval approximation, the joint probability density 
function to find borders in positions x±,X2,--- ,xn in a line of length L is 
given by 



P N (xt,--- ,x N ) = -^P {0 \x 2 -x 1 ) ■ ■ ■p ( - \xn-xn- 1 )p^ ) {x 1 + L-x n ), (22) 
in compact form 

1 N 

P N (x!, ■ ■ ■ , x N ) = — *[[pW(x i+1 - Xi). (23) 



Zn . 
i=i 



where we have considered periodic boundary conditions and we defined xn+i 
X\ + L. The partition function Zn is the normalization constant 



f N 
Z N = / dxi . . . dx N Y\_P { °\xi+i - %i) . (24) 

J x\<X2<---<xm<x\+L j =1 

Comparing ( l23l) with a Boltzmann factor with an inverse temperature (3 = 1, 
we see that, under the independent interval approximation, the statistics of 
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the domain edges is equivalent to a system with N particles which interact 
according to the potential 



N 




i=l 



Thus, it is equivalent to a statistical equilibrium system of particles on a 
circle interacting through a nearest neighbor pair potential. The partition 
function and correlations can be computed by means of the Laplace transform 
as explained in general for this kind of systems in Ref. [IB] . We already 
computed the correlation function g(r) in Eqs. (|20|) and (l2"Tl) . The partition 
function can be obtained as follows. Doing a change of variable 
in Eq. (I24p we have 

Z N = L [ dx 2 ---dx N p(°\x 2 ) m^^-i^lpC^L-^). 

J0<X2<-<x N <L \i=3 J 

(26) 

This is an iV-fold convolution product of p(°>(x). Thus the Laplace transform 
of the partition function is 

J\- lL ^dL=$°\l)) N = fl-le l2 ^eTicf^jY . (27) 

We recall that the circular orthogonal ensemble is equivalent to a system 
of particles on a circle interacting through a log-potential, see Eq. (11 II) . From 
the equations fl25l) and ffTTl) . it is clear the difference between the independent 
domain approximation and the circular orthogonal ensemble, in the first 
system there are only nearest neighbor interactions, while in the second one 
every single particle interact with all the other particles, the pair potential 
is not restricted to nearest neighbors. Additionally, the functional form of 
both potentials is different. 

In figure [5j we compare the gas and spin system correlation functions and 
spacing distributions with the theoretical predictions from the independent 
interval approximation. We notice that this approximation reproduce more 
closely the spacing distributions and the correlation function. The IIA cor- 
relation function has an oscillatory behavior near at r = 1 as it occurs in 
the gas and spin system. Although the predictions are not identical, the IIA 
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gives a good approximation for the correlation function and spacing distri- 
butions with an error of a few percents (5%) deviation from the numerical 
results for the gas and the spin system. 

The fact that the independent interval approximation reproduce much 
better the correlation function and spacing distributions of the gas and the 
spin system than the other approaches considered (CRW, IRW and COE), 
suggest that in the gas and spin systems the domains are not strongly corre- 
lated. 
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Figure 5: Comparison between the gas, spin system and IIA for a lattice 
with N = 1000 sites. 



6 Conclusion 

We studied two out of equilibrium systems, the gas and spin system, which 
present a formation of growing domains. We studied the statistical properties 
of the edges of the domains in the scaling regime and compared them with 
several other models. Since in the scaling regime the statistical properties of 
the edges of the domains are time independent, when properly rescaled, we 
tried to find if there exists an equilibrium statistical system which reproduce 
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these statistics. A first candidate was the circular orthogonal ensemble of 
random matrices (COE), since the scaling size distribution of the domains 
satisfies the Wigner surmise. 

Studying the other spacing distributions p( n \s) and the correlation func- 
tion, we found that only the nearest neighbor distribution (n = 0) coincides 
for the gas system, the spin system, the COE and the models of coalescing 
and interacting random walkers (CRW and IRW) which were introduced in 
Ref. [I] to explain the properties of the gas system. The distribution functions 
p( n \s) for n > 1 and the correlation function are different for each of these 
systems. The nearest neighbor distribution for the coalescing random walk 
is obtained by solving the diffusion equation with an absorbing condition in 
s — [U and the scaling solution turns out to be given by Eq. ([2]). As a 
side note, it is interesting to notice that, up to a normalization, this is the 
wave function of the first excited state of the harmonic oscillator, an analogy 
that has been used in the context of persistence, see Ref. [T9l [20] - Our work 
shows that it is merely a coincidence that this nearest neighbor distribution 
turns out to be the Wigner distribution from the random matrices theory. 

We found a model which provides a better approximation for the gas and 
the spin system domain statistics. This model is the independent interval 
approximation, even though in this approximation the correlation between 
domains is not taken into account as it occurs in interacting and coalesc- 
ing random walk. The independent interval approximation is statistically 
equivalent to an equilibrium system of particles with only nearest neighbor 
interactions. 

The fact that the Wigner surmise, Eq. ([2]), reproduces correctly the near- 
est neighbor spacing distribution of all the different systems we considered 
here, but that the other spacing distributions differ from system to system, 
seems to show that the nearest neighbor distribution has some kind of uni- 
versality while the other distributions do not. However, it can also be in- 
terpreted as an indication that the nearest neighbor distribution does not 
contain enough information about the statistics of the system: several dif- 
ferent systems could share the same nearest neighbor distribution, while the 
finer details and differences between them are contained in the correlation 
function and the other spacing distributions. This is the most important re- 
sult of this paper, because in many cases complex systems are mapped onto 
more simple systems with the only criteria of the nearest neighbor distribu- 
tion similitude. 
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